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Abstract 

We show that, by correctly selecting the probability distribution function p(s) for a 
particle’s distance-to-collision, the nonclassical diffusion equation can be represented 
exactly by the nonclassical linear Boltzmann equation for an infinite homogeneous 
medium. This choice of p(s) preserves the true mean-squared free path of the system, 
which sheds new light on the results obtained in previous work. 


1. Introduction 


A nonclassical linear Boltzmann equation has been recently proposed to address 
particle transport problems in which the particle flux experiences a nonexponential 
attenuation law M- This nonexponential behavior arises in certain inhomogeneous 
random media in which the locations of the scattering centers are spatially correlated, 
such as in a Pebble Bed reactor core PHI El- 

Independent of these developments, a similar kinetic equation has been rigorously 
derived for the periodic Lorentz gas in a series of papers by Golse (cf. [6] ), and by 
Marklof and Strombergsson mi- Related work has also been performed by Grosjean 
[9], considering a generalization of neutron transport that includes arbitrary path- 
length distributions. 

For the case of monoenergetic particle transport with isotropic scattering, the 
nonclassical linear Boltzmann equation is written as 


dip 

ds 


(x, f2, s) + f2 ■ Vq(x, f2, s ) + S t (s)'0(x, f2, s) 


47T 


r /»oo 

c / S t (s / )'0(x, f2', s')ds'dQ' + Q(x) , 

. Jiir JO J 


( 1 . 1 ) 


where ip is the nonclassical angular flux, c is the scattering ratio (probability of 
scattering), and Q(x ) is a source. Here, the total cross section is a function of the 
path length s (distance traveled by the particle since its previous interaction), such 
that the path length distribution 

p(s) = Z t (s)e-f° Et(s ' )ds ' ( 1 . 2 ) 
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does not have to be exponential. If p(s) is exponential, Eq. (1.1) reduces to the 
classical linear Boltzmann equation 


Cl ■ V0(x, Cl) + E t ^(x, Cl) — — f 0(x, Cl')dCl' + 

47r ./Ur 

for the classical angular flux 

poo 

ip(x,Cl) = / ip(x,Cl,s)ds. 


/ , Q( x ) 
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(1.3) 


(1.4) 


It has been shown ra that, by selecting E t (s) in a proper way, Eq. Ql.l| ) can be 
converted to an integral equation for the scalar flux 


0o( x ) = / 0(x, Cl)dCl 

J 47T 


(1.5) 


that is identical to the integral equation that can be constructed for certain diffusion- 


based approximations to Eq. (1.3) in the hierarchy of the SP’n equations HH. 

The work in this paper shows that this is also the case for the nonclassical diffusion 
equation |T| 

/ 9 \ 

C 0o(x) = Q{x), (1.6) 


^ ^ V 2 0o( x ) + 1 


6(s) 


<*> 


which is an asymptotic approximation of Eq. (1.1) when < oo. Here, 


(s) = / sp(s)ds and (s 2 ) = / s 2 p(s)ds. 


(1.7) 


Specifically, we End p(s) and the corresponding E t (s) such that the integral equation 


for Eq. (1.6) is identical to the integral equation for Eq. (1.1). We also show that the 


second moment of p(s) for nonclassical diffusion preserves the true mean-squared free 
path (s 2 ), which gives a new insight on the results obtained in pLOj for the classical 
diffusion approximations. 

The remainder of this paper is organized as follows. In section [2] we convert Eq. 
(1.1) to an integral equation for the scalar flux given in Eq. (1.5). In section [3] we 
convert Eq. (1.6) to an integral equation for the scalar flux, and find the correct 


choice of p(s) such that the integral equation obtained in section [2] is identical to 
this nonclassical diffusion integral equation. We also present a numerical example 
illustrating the differences on p(s) and E t (s) for classical transport, classical diffusion, 
and nonclassical diffusion. The paper concludes with a discussion in section [4j 

2. Integral Equation Formulation 

Let S(x) be given by 


S(x) = c / / E t (s / )0(x, Cl', s')ds'dCt’ + Q(x) 

J 4ir JO 

= c/(x) + Q(x), 


(2.1a) 


2 














where 


and 


/(x) = / S f (s / )0o(x, s')ds' = collision-rate density 

Jo 


(fro (x, s) = / -0(x, f2, s)dQ = nonclassical scalar flux. 

J 47r 


We can now write Eq. (1.1) as the initial value problem 


dxb 

— (x, n,s) + ft ■ V^(x, Q, s) + E t (s)^(x, O, s) = 0, 0 < s, 
as 


^(x,n,o) = 


5(x) 

47T 


(2.1b) 

(2.1c) 

(2.2a) 

(2.2b) 


Following the steps presented in [T] and ra, we: (i) use the method of charac¬ 


teristics to calculate the solution of Eqs. (2.2); (ii) operate on this solution by 
f 47r f 0 T, t (s)(-)dsdQ; and (iii) perform the change of spatial variables from the 3-D 
spherical (O, s) to the 3-D Cartesian x' defined by x' = x - sfl. This yields 


/ (x) = 



MttIx'-xI 2 


(2.3) 


where |x 7 — x|) and S'(x) are given by Eqs. (1.2) and (2.1a), respectively. 


3. Nonclassical Diffusion 


The nonclassical diffusion formulation presented in Eq. (1.6) is an asymptotic 


approximation of Eq. (1.1) in the case of <^s 2 ) < oo. In this formulation, (s') and (s 2 ) 
as given by Eq. ( 1.7[ ) represent the first and second moments of the true path length 
distribution p(s). If p(s) is exponential, then (s ) = 1/S 4 , (s 2 ) = 2/E 2 , and Eq. (1.6) 
reduces to the classical diffusion equation 


—V 2 0 o (x) + (1 - e)£ t 0o(x) = Q(x). 


(3.1) 


For the general case in which p(s) is not assumed to be an exponential, we define 

-l 


S(x) = c(s ) (fro (x) + Q(x) and rewrite Eq. (1.6) as: 

-V 2 0o(x) + A 2 0 o (x) = A 2 (s)5’(x), 

where 
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A = 


(s 2 )- 


(3.2a) 


(3.2b) 


The Green’s function for the operator on the left hand side of Eq. (3.2a) is: 

g—A|x—x'| 


Gdx-x 7 !) = 


47t|x 


(3.3) 
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therefore, we can transform Eq. (3.2a) into an integral equation for 0o( x ) by taking 



(3.4) 


Bearing in mind that (s) represents the mean free path of a particle (i.e. the 
average distance between collisions), the collision-rate density can be written as 
/(x) = (s) 1 0 o ( x ), such that 


/( x ) = 7T = 

< s ) 



\2| _ /| — A|x—x'| 

1 7 | 1 /|2 S(x')dV'. 

47t|x — xy 


This result agrees with Eq. (2.3) iff 


p(s) = \ 2 se Xs = 


6se 


-y/d/<s 2 >s 


b 2 > 


(3.5) 


(3.6) 


It is easy to verify that 


p(s)ds = 


6se 


-a/sT^ 


s /1 >s 


(, 2 > 


-ds = 1, 


(3.7) 


which shows that Eq. (3.6) is a distribution function. Furthermore, E t (s) is given by 

A 2 s 


s,W = 


p(s) 


(3,8) 


jy p(s')ds l 1 + As 

This shows that the nonclassical transport equation reproduces the nonclassical dif¬ 
fusion approximation given by Eq. (1.6) if p(s ) and E t (s) are defined by Eqs. (3.6) 
and (3.8). Moreover, if p(s) is exponential, this results agrees with the p(s) and E t (s) 


obtained for the classical diffusion equation in [[TO] . 
We also point out that 


sp(s)ds = 


s 2 p(s)ds = 


00 6s 2 e - \/ 6 ' /<s2>s 

W 

00 6s 3 e _ \/ 6 / <s2>s 


-ds = 


6 (s 2 ) 


(* 2 > 


ds = (s 2 ). 


(3.9a) 

(3.9b) 


The first moment of p(s) only approximates the mean free path <(s), as can be seen 
in Eq. (3.9a). However, Eq. (3.9b) shows that the true mean-squared free path 
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(]s 2 ) is preserved. This sheds new light on the results presented in [ID], where the 
second moments of the path length distributions obtained for all (classical) diffusion 
approximations give the exact (classical) transport value 2/S 2 . 

Figures [l] and [2] show the functions S t (s) and p(s) for the transport of neutrons 
taking place in the interior of a homogenized 3-D Pebble Bed reactor core system, 
as described in [5j. The system consists of fuel spheres (pebbles) randomly packed 
in a background void with packing fraction 0.5934. The parameters for the material 
of the fuel pebbles are: diameter d — 1; total cross section E t = 1; and scattering 
ratio c = 0.99. As discussed in detail in [5], the classically homogenized system (using 
the Atomic Mix model) has total cross section E* = 0.5934F/ = 0.5934, and its true 
mean-squared free path is numerically calculated to be (s 2 ) = 6.2898. We note that 
nonclassical transport occurs due to the spatial correlations of the fuel pebbles in the 
system, and therefore (s 2 ) ^ 2/E t . 


4. Discussion 


We have shown that the nonclassical diffusion equation for an infinite homogeneous 
medium can be represented exactly by the nonclassical linear Boltzmann equation 
with the correct choice of E f (s) and p(s). We derived explicit expressions for these 
quantities and showed that, while the first moment of the path length distribution 
p(s) only approximates the true mean free path ^s), its second moment preserves the 
true mean-squared free path (s 2 ). This result provides a deeper understanding of the 
results presented in [10] for the second moment of the path-length distributions for 
classical diffusion approximations. 

The work on this paper allows us to construct p(s) and E t (s) that yield the correct 


solution for the nonclassical Boltzmann equation (1.1) in a diffusive system. The only 


parameter necessary for this is the true mean-squared free path (s 2 ). This paves the 
road to the possibility of using this easy-to-obtain p(s) to approximate the solutions 
of the nonclassical Boltzmann equation as the system moves away from the diffusive 
limit. Further work needs to be done to investigate how well such approach would 
perform; this task, however, must be left for future work. 
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Figure 1: Total cross section as a function of s 
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Figure 2: Path length distribution function p(s) 









